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Cosmic microwave background (CMB) polarimeters aspire to measure the faint B-mode signature 
predicted to arise from inflationary gravitational waves. They also have the potential to constrain cosmic 
birefringence, rotation of the polarization of the CMB arising from parity-violating physics, which would 
produce nonzero expectation values for the CMB’s temperature to B-mode correlation ( TB ) and B-mode to 
B-mode correlation (Eli) spectra. However, instrumental systematic effects can also cause these TB and Eli 
correlations to be nonzero. In particular, an overall miscalibration of the polarization orientation of the 
detectors produces TB and EB spectra which are degenerate with isotropic cosmological birefringence, 
while also introducing a small but predictable bias on the BB spectrum. We find that BicepI three-year 
spectra, which use our standard calibration of detector polarization angles from a dielectric sheet, are 
consistent with a polarization rotation of a = — 2.77° ± 0.86°(statistical) ± 1.3°(systematic). We have 
revised the estimate of systematic error on the polarization rotation angle from the two-year analysis by 
comparing multiple calibration methods. We also account for the (negligible) impact of measured beam 
systematic effects. We investigate the polarization rotation for the BicepI 100 GHz and 150 GHz bands 
separately to investigate theoretical models that produce frequency-dependent cosmic birefringence. We 
find no evidence in the data supporting either of these models or Faraday rotation of the CMB polarization 
by the Milky Way galaxy’s magnetic field. If we assume that there is no cosmic birefringence, we can use 
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the TB and EB spectra to calibrate detector polarization orientations, thus reducing bias of the cosmological 
R-mode spectrum from leaked E-modes due to possible polarization orientation miscalibration. After 
applying this “self-calibration” process, we find that the upper limit on the tensor-to-scalar ratio decreases 
slightly, from r < 0.70 to r < 0.65 at 95% confidence. 

DOI: 10. 1 103/PhysRevD.89. 062006 PACS numbers: 98.70. Vc 


I. INTRODUCTION 

The cosmic microwave background (CMB) is a powerful 
cosmological probe; recombination physics, structure for- 
mation, and the cosmological reionization history represent 
only a small subset of the phenomena probed by its 
temperature and polarization anisotropy. In addition, sev- 
eral aspects of fundamental physics can be constrained by 
CMB observations, the most familiar of which are infla- 
tionary physics revealed via the imprint of primordial 
gravitational waves in the polarization of the CMB and 
the masses of neutrinos which can be probed via gravita- 
tional lensing by dark matter. These phenomena create 
E-mode polarization at the sub-// A" level. 

Cosmological information can be extracted from the 
CMB’s power spectra. Out of the six possible pairings of 
the temperature anisotropy T and polarization E- and 
E-modes, only four have nonvanishing expectation values 
in the standard (ACDM) cosmological paradigm. The 
expectation values of the TB and EB cross correlations 
vanish in the standard cosmological model due to parity 
symmetry but may assume nonvanishing values in the 
presence of systematics, astrophysical foregrounds, or, more 
interestingly, parity-violating departures from the standard 
models of electromagnetism and gravity. Any mechanism 
capable of converting E- to E-mode polarization necessarily 
leaks the temperature to E-mode correlations (TE) and 
E-mode auto spectrum (EE) to TB and EE, respectively. 

A detection of EE and EE correlations of cosmological 
origin could undermine the fundamental assumptions of 
parity symmetry and Lorentz invariance by showing that 
our Universe possesses a small degree of chirality. This 
phenomenon can be best revealed by CMB polarization 
where minuscule effects can accrue to observable levels over 
the 13.8 Gyrs since CMB photons last scattered from the 
primordial plasma. This preferred chirality can be induced 
by the coupling of a pseudoscalar field to either Chem- 
Simons-type terms in the electromagnetic interaction [1-3] 


or the Chem-Pontrayagin term in the case of gravitational 
interactions [4-6]. This work constrains the parameters in a 
scale-independent cosmological birefringence model as well 
as investigating frequency-dependent scale-independent 
models. Current best constraints (not including this work) 
on scale-independent cosmological birefringence from CMB 
experiments are shown in Table I. Though constraints on 
scale-dependent birefringence models have been reported 
with WMAP data [7-9], we do not provide such constraints 
in this work. A 3/7 detection of cosmic birefringence was 
reported from combined WMAP, BOOMERanG, and 
BicepI two-year results (while explicitly excluding QUaD 
data) in [10]. This work was later updated to include the 
impact of systematic effects at the levels reported by the 
three experiments and the significance reduced to 2.2/r [11]. 

BicepI has set the most stringent constraints on the 
CMB’s E-mode power spectrum [16,17] in the multipole 
range 30 < € < 300. BicepI also measured the TB and EE 
power spectra in this range [16,17]. These TB and EE 
modes are extremely sensitive probes of departures from 
the standard cosmological model. In this work, we analyze 
the full BicepI three-year spectra [17] for evidence of 
polarization rotation, considering systematic uncertainties 
including our primary and alternate polarization calibra- 
tions, and exploring constraints on cosmological birefrin- 
gence. We then use this polarization angle to self calibrate 
detector polarization orientations and calculate the tensor- 
to-scalar ratio from the self-calibrated EE spectrum. 

The outline of the paper is as follows: Sec. II contains a 
review of polarization rotation and how it affects the 
observed CMB power spectra. The data sets and analysis 
procedure are described in Sec. III. Results and consistency 
checks are presented in Sec. IV. The impact of instrumental 
systematics is discussed in Sec. V. Consistency of the data 
with different birefringence models is in Sec. VI. Application 
of self calibration and its effect on the tensor-to-scalar ratio, r, 
are in Sec. VII, and we discuss our results in Sec. VIII. 


TABLE I. Previous rotation angle constraints from CMB experiments, following [12]. Systematic uncertainties 
are shown in parentheses, where provided. 


Experiment 

Frequency (GHz) 

f range 

a (degrees) 

WMAP7 [13] 

41+61+94 

2-800 

-1.1 ± 1.4(±1.5) 

BOOM03 [14] 

143 

150-1000 

-4.3 ±4.1 

QUaD [15] 

100 

200-2000 

—1.89 ± 2.24(±0.5) 

QUaD [15] 

150 

200-2000 

+0.83 ± 0.94(±0.5) 
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II. POLARIZATION ROTATION OF THE CMB 
POWER SPECTRA 

The CMB can be described by the statistical properties of 
its temperature and polarization. E- and ri-rn ode polariza- 
tion can be formed from linear combinations of the Stokes 
Q and JJ parameters. Maps of the temperature, T, and 
Stokes parameters, Q and U, are expanded in scalar and 
spin ±2 spherical harmonics [18,19] to obtain 

T {n) = Y^a T t m Y fm {h) 

(Q ± iU){h) = 5> ±2/m±2 Wn), (1) 

f,m 

where the E- and B-modes of polarization have expansion 
coefficients a E m and a B m which can be expressed in terms 
of the spin ±2 coefficients 

a±2,fm = a 'L ± ia 'L- ( 2 ) 

The spherical harmonic coefficients, a^ m , are characterized 
by their statistical properties, 

<4J = o 

( a tm a fm') = C X /'5^5mm', (3) 

where X and X' are either T. E or B. Here, (a) stands for the 
ensemble average. The polarization modes, E and B. are 
pure parity states (even and odd, respectively) and thus the 
correlation over the full sky of the /I -mode with either the 
temperature or E-mode polarization vanishes [18,19]. 
However, if the polarization of the CMB is rotated, there 
will be a mixing between E and B. subsequently inducing 
TB and EB power spectra (Fig. 1): 
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r'TT _ r TT 

C'J E = C EE cos(2a) 

C' EE = C EE cos 2 (2a) + C BB sin 2 (2a) 

C IBB = Cf £ sin 2 (2cr) + C| B cos 2 (2a) 

C'J B = C EE sin(2a) 

C' EB = i ( Cf E - Cf B ) sin(4 a ) . (4) 

No assumption has been made here as to the source of this 
rotation, namely whether or not it is cosmological. In the 
literature, a is identified with the birefringence rotation 
angle (see [13,15]), though here it is used to denote 
polarization rotation of any origin. 

III. DATA AND ANALYSIS METHODOLOGY 

BicepI observed for three years at the South Pole in 
three frequency bands: 100, 150 and 220 GHz, and released 
two-year results from 100 and 150 GHz frequency- 
combined spectra in [16] and three-year frequency- 
combined spectra in [17]. Results from the BicepI 100, 
150, and 220 GHz observations of the galactic plane are in 
[20] and from Faraday rotation modulators in [21]. 

We employ maximum-likelihood estimation for deter- 
mining the best-fit polarization rotation angles of the power 
spectra following Eq. (4). We use two methods to construct 
the likelihoods, a Gaussian bandpower likelihood approxi- 
mation and the Hamimeche-Lewis likelihood construc- 
tion [22]. 

A. Data sets 

We calculate rotation angles from the three-year fre- 
quency combined “all-spectra” estimator, where all-spectra 
is defined as TE + EE + BB + TB + EB. We can break 
this down by frequency and by spectral estimator for 
consistency checks. From this, we get four frequency 



FIG. 1 (color online). Standard ACDM power spectra after applying polarization rotation of —3° (blue) to +3° (red), in 0.5° steps, for 
TB (left) and EB (right). 
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subsets consisting of the two frequency autospectra: 
100 GHz autospectra (denoted “100”) and 150 GHz 
autospectra (denoted “150”), and the two frequency cross 
spectra: 100 GHz cross correlated with 150 GHz (denoted 
“cross”) and 150 GHz cross correlated with 100 GHz 
(denoted “alt-cross”). Note that although the EE and BB 
spectra are identical for the cross and alt-cross data sets, the 
TB and EB spectra are not, e.g., T l00 B 150 ^ j^sogioo. pj ots 
of the T B and EB spectra for these four frequency subsets 
are in Fig. 2. 

In addition, we have four spectral combinations to 
constrain a: the TB and EB modes as well as the 
combination of TB + EB, and all-spectra: TE + EE + 
BB + TB + EB since polarization rotation also affects 
TE. EE, and BB: however, from Eq. (4), we can see that 
for small a the rotated TE, EE, and BB deviate from the 
unrotated spectra by order a 2 and thus their constraining 
power for a is much weaker than the T B and EB spectra, 
which are linear in a. In addition, since they are quadratic in 
a, the sign of a cannot be directly determined. TB or EB 
break this sign degeneracy. These are not independent 
estimators but are useful as any unexpected discrepancies 
can be used to test the validity of the analysis. 



0 50 100 150 200 250 300 350 


I 

(a) TB 



FIG. 2 (color online). BicepI TB and EB power spectra for all 
frequency combinations: 100 GHz autospectra (red), 150 GHz 
autospectra (blue), 100 x 150 GHz cross-spectra (green), 150 x 
100 GHz “alt-cross" spectra (cyan), and frequency combined 
100 + 150 GHz spectra (black). The points have been displaced 
in € for clarity. 
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B. Likelihood analysis 

We employ two likelihood constructions for this analy- 
sis: a Gaussian bandpower likelihood approximation and 
the more accurate Hamimeche-Lewis (HL) likelihood 
approximation [22]. The two likelihood constructions 
produce similar results, although we use the HL method 
for the final results since it more accurately treats cross- 
spectra covariances. We test both likelihood constructions 
for any biases and, in simulations, we find they accurately 
recover known input rotation angles. 

For both methods, we calculate/ 2 = — 21n£, where/ 2 
is defined in Eq. (6), below. We found the rotation angle 
that maximized the likelihood, and constructed la error 
bars by finding the minimum-width 68% credible interval, 
assuming a uniform prior on a, for both likelihood 
constructions. 

1. Gaussian bandpower likelihood approximation 

This method was chosen due to its computational 
efficiency for isolating individual spectral estimators with- 
out including corresponding autospectra. Here, the differ- 
ence between the observed spectra and theory spectra 
including rotation, 

A f(a) = Pf -X>f(a), (5) 

is computed as a function of rotation angle, a, for each 
BicepI multipole bin, where BicepI reports nine bins of 
uniform width Af = 35, with the first bin spanning 20 < 
€ < 55 and the ninth bin spanning 300 < € < 335. Here, 
V? is the measured BicepI XY power spectrum and 
Vf/ (a) is the theoretical rotated bandpower for XY for a 
given a. We use T>y to denote binned estimates of 
vy = €(€ + \)cy/2n. Here, XY = TB or EB for each 
frequency combination. The / 2 statistic is then constructed 
using 

Xxy(<*) = '^ A b Y ( a ) M bl' A b' Y ( a ), (6) 

bV 

where M. bb ' is the covariance between multipole bins b and 
b' , as described in [17]. 

2. Hamimeche-Lewis method 

The Hamimeche-Lewis method is the bandpower like- 
lihood approximation used in [17]. As before, the / 2 
statistic is constructed as in Eq. (6) but following the 
procedure outlined in [17]. One crucial difference between 
this method and the Gaussian bandpower likelihood 
approximation is that XY includes all combinations of 
the spectra X and Y. For example, for EB, this method does 
not calculate the / 2 for EB but actually the / 2 which 
includes EB + EE + BB — the comparison of the measured 
spectra to theoretical rotated spectra for EB, EE and BB 
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simultaneously. To calculate the/ 2 statistic for any “pure” 
spectral combination using this method, we calculate the/ 2 
of the full spectral combination and subtract off the 
other spectral combinations. For example, for EB, /| B = 

Xee+bb+eb Xee Xbb- 

IV. ROTATION ANGLE RESULTS 

The rotation angle, a , was calculated using the HL 
method from the standard Bicep 1 three-year frequency 
combined spectra. The best-fit rotation angle is a = 
—2.77° ± 0.86°, where the quoted uncertainty is purely 
statistical. These spectra use our standard calibration of 
detector polarization angles from a dielectric sheet; sys- 
tematic uncertainty on this calibration is discussed below in 
Sec. V. Figure 3 plots the peak-normalized F1L likelihood 
and Fig. 4 shows the best-fit rotation angle spectra plotted 
compared to the Bicep 1 three-year data and 499 simulated 
ACDM realizations (i.e., with a = 0). 

A. Consistency between analysis methods 

To check that the rotation angle is not dependent on the 
analysis method, polarization rotation angles derived from 
the two analysis methods were compared. In Fig. 5, the 
likelihoods calculated for T B, EB and TB + EB for the 
frequency-combined spectra for both the HL and Gaussian 
bandpower likelihood approximations are overplotted. For 
all three available spectral estimators, the analysis methods 
agree to within 0.32c, 0.30c and 0.18c for TB, EB and 
TB + EB, respectively. 



FIG. 4 (color online). Frequency-combined three-year BicepI 
spectra (black points) shown with the theoretical rotated spectra 
from the best-fit all-spectra rotation angle, a = —2.77° ± 0.86° 
(red solid), the lc confidence limits (red dashed), and the 499 
simulation realizations (gray). All simulation realizations 
assume a = 0. 


B. Consistency between frequencies 

For consistency, the different frequency combinations 
were checked to determine if they have similar rotation 
angles. Table II shows the calculated rotation angles from 



a - rotation angle (deg) 

FIG. 3. The peak-normalized Hamimeche-Lewis likelihood for 
the all-spectra (“ TEB ”) rotation angle. The maximum likelihood 
value is — 2.77 ° and the 68% confidence limits are ±0.86° from 
the peak value, corresponding to 3.22(7 statistical significance. 


each frequency data set and for all four spectral estimators. 
Figure 6 shows the HL likelihoods for the all-spectra (T EB) 
rotation angles for each data set. 

C. Consistency with Planck temperature taps 

The TB spectrum estimate of the polarization rotation 
explicitly depends on the BicepI measurement of temper- 
ature. To check for systematics in the TB power spectrum, 
we replace BicepI maps with Planck temperature maps 
[23] for both 100 and 150 GHz and find the recovered 
angles agree to within 0.2(7. 

V. IMPACT OF INSTRUMENTAL 
SYSTEMATIC EFFECTS 

BicepI was the first experiment designed specifically to 
measure the /(-mode power spectrum in order to constrain 
the inflationary cosmological model [24]. Accordingly, the 
analysis of instrumental systematics focused on potential 
bias of the BB spectrum and the tensor-to-scalar ratio with a 
benchmark of r = 0.1 [25]. Here, we extend the analysis to 
include the impact of measured systematics on the T B and 
EB power spectra for the three-year data set. 


062006-5 


J. P. KAUFMAN 



FIG. 5 (color online). Comparison of the two likelihood 
methods employed for the rotation angle calculation for the 
frequency-combined data. Likelihoods computed from the 
Gaussian bandpower likelihood approximation are shown as 
dashed lines and the Hamimeche-Lewis method likelihoods 
are shown as solid lines. The TB likelihood is in blue, EB is 
in green, and the TB + EB likelihood is cyan. The likelihoods 
have been normalized by their integral over a. 

A. Polarization angle calibration 

An error in the detector polarization angles used for map 
making is the only systematic which is completely degen- 
erate with a rotation due to isotropic cosmic birefringence, 
and the only systematic capable of producing self- 
consistent TB and EB power spectra [26]. This calibration 
requirement is much more stringent when attempting to 
measure a than for r. 

Calibrating detector angles for CMB polarimeters is very 
challenging. Some commonly employed methods include 
man-made calibrators, such as polarizing dielectric sheets 
[24,25] or polarization-selecting wire grids [27,28], and 
observations of polarized astronomical sources [29,30]. 
Man-made polarization calibration sources suffer from a 
host of challenges: they are often situated in the near field 
of the telescope, they can be unstable over long time scales, 
and they can be cumbersome to implement and align. 
Astronomical sources are not visible from all observatories 
and even the best characterized sources have orientations 
measured to an accuracy of only 0.5° [31]. In addition, the 
brightness of both astronomical and man-made calibration 
sources can overload the detectors, forcing them into a 
nonlinear response regime [32]. 

BicepI employed several hardware calibrators to mea- 
sure detector polarization angles. The primary calibration 
comes from a dielectric sheet calibrator (DSC), described in 
detail in [25], but additional calibrations were made using 
sources with polarizing wire grids in the near and far field. 
The BicepI beam size and observatory location prevented 
polarization calibration using astronomical sources. 

The polarization angle measurement from the DSC was 
performed the most frequently and is the best studied, 
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TABLE II. Maximum likelihood value and la error for a. All 
numbers are in degrees. 


Dataset 

TB only 

EB only 

TB + EB 

All-spectra 

100 GHz 

— 1 7Q+ 3 - 18 
' — 3.14 

q co+2.38 
J,JJ -2.26 

9 99+2.O6 
+•+ ' -1.98 

9 99+2.O6 
+•+'-2.02 

150 GHz 

—4 37+1-92 
7 — 1.78 

—2 95+ 1 - 20 
+*- 7 -1.18 

0 1 0 + 1.14 
— 1 .12 

—2 Q1 + 106 
+*''1 — 1.04 

Cross 

_3 Q3+1.84 

9 cc+1.68 
+* JJ -1.60 

— 2.83+ 1 2 2 4 8 

—2 67 +L20 
+- 0/ — 1.18 

Alt-cross 

—2 71 + 3 - 52 
+- /1 -3.74 

9 9 c +2.26 
2.20 

—3.45+|.f4 

-3 15 +1 - 96 
•~ , - 1J - 2 .oo 

Comb 

—3 47+1-66 
-1.56 

_3 05+1.00 

J.UJ-0 96 

-2 9Q +0 - 94 
+ , ' / ~-0.92 

9 99 + 0.86 

+•' '- 0.86 


which is why it was chosen for results in [16,17], as well as 
this work. Repeated measurements during each observing 
season produced polarization angles that agree with an 
root-mean-square (rms) error of 0.1°. However, measure- 
ments taken before and after focal plane servicing between 
the 2006 and 2007 observing seasons show an unexplained 
rotation of 1° in the polarization angles. There is also some 
uncertainty in translating the results of the DSC measure- 
ment to parameters appropriate for CMB analysis. The 
details of the polarized signal from the dielectric sheet 
depend on the near-field response of each detector, which is 
not well characterized. 

We also consider two alternate calibrations for the 
detector polarization angles, which were both described 
in [25] as methods to measure the cross-polar response of 
the detectors. The first is a modulated broadband noise 
source, broadcasting via a rectangular feedhorn located 
behind a polarizing wire grid. The source is located on a 
mast at a range of 200 meters. We measure the detector 
response as a function of angle by scanning over the source 
with 18 different detector orientations. The advantage of 
this method is that the source is in the far field of the 
telescope. A challenge is that the observations require the 
use of a flat mirror, complicating the pointing model. In 
addition, it takes a significant amount of time to perform 
scans at all 18 orientations, which makes it more difficult to 
maintain stable source brightness. For Bicep 2, we have 
invested significant effort in improving polarization ori- 
entation calibrations with the far-field broadband noise 
source, both by developing a high-precision rotating 
polarized source and improving the pointing model used 
for calibration analysis. These improvements were moti- 
vated by BicepI experience with systematic uncertainties 
on both the DSC and broadband source calibrations. 

Another polarization angle calibrator consists of a wire 
grid covering a small aperture that is chopped between an 
ambient temperature absorber and cold sky. For calibra- 
tions, this source is installed in the near field of the 
telescope and the wire grid is rotated to measure detector 
polarization angles. The interpretation of results from this 
source has significant uncertainty because the small aper- 
ture probes only a small fraction of the detector near field 
response, yet the results are extrapolated to the full beam 
response. 
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FIG. 6 (color online). Comparison of integral-normalized 
Hamimeche-Lewis method likelihoods of all-spectra (TEB) 
rotation angles for 100 GHz autospectra (red), 150 GHz auto- 
spectra (blue), 100 x 150 GHz cross spectra (green), 150 x 
100 GHz “alt-cross” spectra (cyan), and frequency-combined 
spectra (black). 

Table III lists the values of a measured from maps made 
using each of the polarization angle calibrations. Also 
shown is the result obtained if we simply assume that the 
detector polarization angles are as designed. These derived 
a values are qualitatively consistent with the average 
difference in the detector polarization angles between 
any two calibration methods, though the details depend 
on how each detector is weighted in the three-year maps. 
Besides the global rotations between each calibration 
method, which contribute to the variation in a, the per- 
detector polarization angles show scatter of 0.6-0. 9° 
between methods, much larger than the 0.1° consistency 
seen from repeated measurements using the DSC. Despite 
this scatter, we can observe significant structure in the 
pattern of polarization angles from detector to detector, 
which is not present in the as-designed angles. 

From consideration of the 1.14° difference between 
alpha as derived from the DSC calibration and the mean 
of the three alternate calibrations, which have very different 
sources of systematic uncertainty, as well as the 1° shift 
observed in the DSC calibration results between observing 
seasons, we assign a calibration uncertainty of 1.3° on the 
overall orientation from the DSC calibration. We believe 
this upward revision of the 0.7° uncertainty quoted for this 
same calibration in [25] is warranted by the tension with the 
alternate calibrations. While this systematic error is larger 
than the 0.86° statistical error on a, we stress the fact that 
the polarization angle calibration is quite a bit better than 
what is needed to meet the r = 0. 1 benchmark for the 
primary BicepI science goal. 

In Sec. VII, we adopt a different approach and self 
calibrate the polarization orientations by rotating the 
polarization maps to minimize a. Note that the calibration 
uncertainty on a applies only when we use the DSC 
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TABLE III. Polarization rotation angles derived using different 
detector polarization angle calibrations: the dielectric sheet 
calibrator (DSC), the far-field wire grid broadband noise source, 
the near-field wire grid aperture source, and assuming the 
polarization angles are as designed. For all cases, the statistical 
uncertainty on a is 0.86°. 


Calibration method 

Near/far field 

a (degrees) 

DSC 

near 

-2.77 

Wire grid broadband source 

far 

-1.71 

Wire grid aperture source 

near 

-1.91 

As designed 

— 

-1.27 


calibrated maps and attempt to measure astrophysical 
polarization rotation. To judge how well the self-calibration 
procedure can debias the B-mode map, only the statistical 
error is relevant. 

B. Differential beam effects 

Differential beam mismatches potentially mix C-ruodes 
and B-modes or leak intensity to either E- or B-mode 
polarization. Here, we investigate the impact of differential 
beam size, differential relative gain, differential pointing, 
and differential ellipticity on the derived rotation angle. 

Beam systematics affect the EB spectra in a different 
way than TB spectra [33,34]. As a result, the scale 
dependence of the beam systematic polarization will imply 
a different effective rotation angle in the TB spectrum 
versus the EB spectrum, for a fixed range. 

The BicepI beams were measured in the lab prior to 
deployment using a source in the far field (50 meters from 
the aperture) and each observing season during summer 
calibration testing. Beam maps were fit to a two-dimensional 
elliptical Gaussian model which included a beam location, 
width, ellipticity, and orientation of the major axis of the 
ellipse with respect to the polarization axes. 

1. Differential beam size 

Though the differential beam size effect can leak temper- 
ature to polarization, due to circular symmetry it will not 
break the parity of the underlying sky and thus cannot 
generate the parity-odd TB and EB modes [33]. 

2. Differential relative gain 

As with differential beam size, circular symmetry is 
preserved by differential relative gain and thus there is no 
breaking of the parity of the sky which would generate TB 
and EB modes [33]. We ran differential gain simulations 
using observed values and random values drawn from a 
Gaussian distribution with an rms of 1%. None of the 
simulations showed polarization rotation greater in magni- 
tude than 0.25°. 

A significant difference between the BicepI two-year 
results reported in [16] and the BicepI three-year power 
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spectra is that the three-year spectra undergo relative gain 
deprojection which reduces 5-mode contamination due to 
this systematic to negligible levels [17]. 

3. Differential pointing 

The effect of differential pointing is analytically calcu- 
lated using the measured magnitude and direction of beam 
offsets with the expected amount of false BB power scaling 
as the square of the magnitude of the differential pointing, 
following the construction in [33]. The upper limit on 
differential pointing error was estimated to be < 1.3% of the 
beam size. While this was found to be the dominant beam 
systematic effect for the BicepI limit on r [17], it is clear 
from Fig. 7 that differential pointing does not induce TB or 
EB, and has a negligible effect on the polarization rotation 
angle estimation. This was calculated for the worst-case 
scanning strategy and therefore provides very conservative 
bounds on the TB and EB produced. 

4. Differential ellipticity 

Differential ellipticity values were derived by fitting each 
measured beam in a detector pair for ellipticity and then 
differencing those values for the two detectors in a pair. The 
fits were generally not repeatable when the telescope was 
rotated about its boresight angle, so only upper limits on 
differential ellipticity are quoted. The BicepI estimated 
differential ellipticity is <0.2%. 

As with differential pointing, we calculate the TB and 
EB following the construction in [33]. As before, this is for 
the worst-case scenario, where the major axes of the 
ellipticities are separated by 45°. From Fig. 7, it is clear 
that differential ellipticity can generate TB power which 
has a different spectral shape than that produced by 


polarization rotation. In addition, the TB spectrum is 
inconsistent with the polarization rotation of EB. 

C. Experimental consistency checks 

To probe the susceptibility of BicepI data to systematic 
effects irrespective of origin, [16] created six null tests or 
“jackknife” spectra that were used as consistency tests. 
These tests involve splitting the data in two halves and 
differencing them. The two halves are chosen to illuminate 
systematics since signals which are common to both data 
sets will cancel, and the resultant jackknife will either be 
consistent with noise or indicate contamination. The jack- 
knife splits are by boresight rotation, scan direction, 
observing time, detector alignment, elevation coverage, 
and frequency, as described in [16,17]. 

In the power spectrum analyses [16,17], the jackknife 
maps were obtained by differencing the maps for each half 
whereas here we ran each jackknife half through the full 
analysis pipeline to produce power spectra. Unlike in the 
previous analyses, we did not look for consistency with 
zero but self consistency between jackknife halves. We fit 
rotation angles for each jackknife half and difference the 
resultant best-fit angles to form A a. Only the frequency- 
combined spectra are used to improve the constraining 
power of these tests. We then calculate the probability to 
exceed (PTE) the observed A a value by chance, given 
measurement uncertainties. These results are shown in 
Fig. 8. If there is an instrumental systematic contribution 
to a detection of a that is tested using these null tests, then 
excess PTE values near zero will arise. The jackknife 
halves are considered consistent if they meet the following 
three criteria: (1) fewer than 5% of the jackknives have PTE 
values smaller than 5%, (2) none of the PTE values are 
excessively small (defined as c 1%), and (3) the PTE value 
from all jackknives are consistent with a uniform 



FIG. 7 (color online). A plot of the effects of differential ellipticity (red solid line) and differential pointing (red dashed line) on the TB 
(left) and EB (right) power spectra. The gray line shows the power spectra with a = —2.11°. The black points are the frequency- 
combined three-year band powers. The differential ellipticity and differential pointing are at the level of 0.2% and 1.3%, respectively, 
corresponding to the upper limits reported in [25]. In both cases, the systematic curves correspond to “worst-case” scenarios: the major 
axes of the ellipticities are separated by 45° and differential pointing assumes a poorly chosen scan strategy [34]. 
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PTE 


FIG. 8 (color online). A histogram of the probability to exceed 
values for the measurement of A a in each of the jackknife spectra. 
There are 20 spectral combinations (which are not all indepen- 
dent): TB, EB, TB + EB, and all-spectra estimators for the 
frequency-combined spectra for each of the five jackknife tests. 


distribution between zero and one. Given the consistency of 
the jackknife PTEs with these criteria, systematics probed 
by these jackknifes are not the source of the observed 
polarization rotation angle. 

VI. CONSTRAINTS ON FREQUENCY- 
DEPENDENT COSMOLOGICAL 
BIREFRINGENCE 

This paper has focused on the assumption that polari- 
zation rotation is independent of electromagnetic fre- 
quency. However, several models feature polarization 
rotation that predicts a manifestly frequency-dependent 
rotation angle. One such birefringence model has been 
proposed by Contaldi, Dowker, and Philpott in [35], 
hereafter called the “CDP” model. 

Another effect which could cause frequency-dependent 
polarization rotation would be Faraday rotation of CMB 
polarization due to the Milky Way’s magnetic field. 

A. Contaldi-Dowker-Philpott model 

In the CDP model, there are two electromagnetic 
frequency-dependent parameters (pi and /) leading to the 
following power spectra: 

r'TT _ r'TT 

C’J B = e^C EE sin(2/) 

C' EB = l -e~ 2 ^C EE - C BB ) sin( 4 /) 

C’J E = e -»C T / cos(2/) 

C’ EE = e- 2f, C E / cos 2 (2/) + e~ 2lJ C BB sin 2 (2%) 

C ,BB = e ~^Cf E sin 2 (2/) + e~ 2 ^C B / cos 2 (2j), (7) 
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TABLE IV. Maximum likelihood values for the damping 
parameter, pi, and the rotation angle, 7 , along with their 1 <t error 
bars for the CDP model. 


Frequency (GHz) 


7 (degrees) 

100 

-0 oi7 +0 073 

U.Ul / —0.076 

9 9 ^ + 2.02 
2.02 

150 

—0 029 +0 042 
u.uz ^_ 0043 

— ? Q 1 + 102 
z ' yi -1.02 


where e/x ~ 1 /g and v is the electromagnetic frequency 
(i.e., 100 and 150 GHz). Here, 7 is a frequency-dependent 
rotation angle, and pi characterizes the frequency-dependent 
damping parameter. The frequency-independent spectra are 
obtained in the limit pi -* 0, with 7 identified with a. 

As is evident from Eq. (7), to constrain the frequency- 
dependent CDP model, the TE, EE and B B spectra must 
also be included in the analysis in order to break the 
degeneracy between 7 and pi. 

The results for the damping parameter and rotation angle 
7 are presented in Table IV and Fig. 9. The inferred pi is 
consistent with zero and pi/% is not inversely proportional to 
v, thus there is no compelling evidence for frequency 
dependent birefringence in the CDP picture. 

B. Faraday rotation of galactic magnetic field 

Faraday rotation due to the Milky Way’s magnetic field 
could produce frequency-dependent polarization rotation 
proportional to the inverse of the frequency squared [36]. 
Scaling the 150 GHz all-spectra estimator a by 
(100/ 150) -2 , we would expect to see polarization rotation 
of the 100 GHz spectra consistent with an a = — 6.55lf 39 
degrees, whereas our 100 GHz all-spectra estimator results 
in a polarization rotation of a = — 2.27// f ( p , corresponding 
to a 1.37(7 discrepancy. While there is some tension 
between the predicted a and the measured a, Faraday 
rotation cannot be ruled out as the cause of the rotation. 





FIG. 9 (color online). Best-fit pi and 7 values for the CDP model 
for 100 GHz (red point) and 150 GHz (blue point) from the all- 
spectra estimator, along with their 68 % confidence interval 
contours (red and blue contours, respectively). 
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VII. SELF-CALIBRATED UPPER LIMIT ON 
TENSOR-TO-SCALAR RATIO 

If the polarization rotation is systematic in nature, the 
derived rotation angle can be used to calibrate the detector 
polarization orientations [32]. The three-year all-spectra 
rotation angles were added to the polarization orientations 
treating the frequency bands as independent, i.e., only the 
100 GHz (150 GHz) derived rotation angle was added to 
the 100 GHz (150 GHz) detectors. These self-calibrated 
polarization orientation angles were propagated through the 
power spectrum analysis pipeline [17]. The self-calibrated 
power spectra were analyzed for residual polarization 
rotation which yielded a rotation angle a = +0.01° ± 
0.86° from the frequency-combined all-spectra estimator, 
consistent with zero, as expected. 

Any polarization rotation, regardless of cosmic or 
systematic origin, will positively bias r since E-mode 
power will be leaked into the E-mode spectrum 
[Eq. (4)]. There is also a reduction to the E-mode power 
spectrum due to E-modes leaking to E-modes; however, 
since the E-modes are significantly larger than the 
E-modes, the net result is a positive bias on the E-mode 
power spectrum. From the self-calibrated three-year power 
spectra, following the procedure in [17], we find the upper 
limit on the tensor- to-sc alar ratio reduces from r < 0.70 to 
r < 0.65 at 95% confidence. From simulations, we find 
that the bias on r from self calibration with no underlying 
polarization rotation is less that 0.01. 

VIII. CONCLUSION 

The Bicep 1 three-year data, when analyzed using 
detector polarization orientations from our standard dielec- 
tric sheet calibrator, show nonvanishing T B and EE spectra 
consistent with an overall polarization rotation of —2.11° ± 
0.86° at 3.22 a significance. The significance for nonzero 
rotation of astrophysical origin is only 1.78n, given the 1.3° 
systematic uncertainty on our orientation calibration which 
adds in quadrature. This result passes experimental con- 
sistency tests which probe for systematic differences of 
polarization rotation in various subsets of data. We rule out 
beam systematics as significant, and identify polarization 
orientation miscalibration as the primary concern among 
instrumental systematics. Isotropic cosmic birefringence 
cannot be excluded, though it is degenerate with a polari- 
zation miscalibration. The data show no compelling 
evidence for frequency-dependent isotropic cosmic bire- 
fringence models. An alternate use of the measurements 
described here is to self calibrate the detector polarization 
orientations, at the expense of losing constraining power on 
isotropic cosmological birefringence [32]. Self calibrating 
the BicepI three-year data reduces the upper limit on the 
tensor-to-scalar ratio from r < 0.70 to r < 0.65 at 95% 
confidence. 


Future CMB polarimeters with improved polarization 
calibration methods will be needed to break the degeneracy 
between polarization rotation and detector polarization 
orientation uncertainty. In addition to the CMB, compli- 
mentary astronomical probes such as the polarization 
orientation of radio galaxies and quasars [37,38] can help 
constrain cosmological birefringence. However, these 
objects can only constrain cosmic birefringence over a 
limited range of redshifts and only along particular lines of 
sight, whereas CMB polarization can be used to constrain 
cosmic birefringence over the entire sky and is sensitive to 
effects accrued over the history of the entire Universe. 
Polarization angles calibrated with current man-made or 
astronomical sources are accurate enough for current 
generation E-mode measurements, but are insufficiently 
characterized for cosmic birefringence searches. Based on 
BicepI experiences with systematic uncertainties on 
polarization orientation calibration reported in this paper, 
improved far-field calibrators have been developed for 
Bicep2 and other future experiments. The revolutionary 
discovery potential of a detection of cosmic birefringence 
motivates the development of more accurate hardware 
calibrators and further investigation of astronomical 
sources to achieve a precision of <sc0.5°. Ultimately, a 
combination of precisely understood man-made and astro- 
nomical sources will allow for powerful constraints on 
parity violation which will come concomitantly with 
bounds on the physics of inflation. 
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